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Summary 

A new, structurally compliant rocket engine combustion 
chamber design has been validated through analysis and ex- 
periment. Subscale, tubular channel chambers have been 
cyclically tested and analytically evaluated. Cyclic lives were 
determined to have a potential for 1000 percent increase over 
those of rectangular channel designs, the current state of the art. 
Greater structural compliance in the circumferential direction 
gave rise to lower thermal strains during hot firing, resulting in 
lower thermal strain ratcheting and longer predicted fatigue 
lives. Thermal, structural, and durability analyses of the 
combustion chamber design, involving cyclic temperatures, 
strains, and low-cycle fatigue lives, have corroborated the 
experimental observations. 


Introduction 

Rocket combustion chambers, such as the one used in the 
space shuttle main engine, operate at extremely high heat flux 
(100 Btu/in. 2 -sec at the throat). The high heat flux results from 
high-temperature (3315 °C; 6000 °F), high-pressure (20 MPa; 
3000 psi) combustion gases. Active cooling of the chamber is 
mandatory to keep the hot-gas-side wall material from melting. 
The choice of materials for constructing the hot-gas-side wall 
is limited. Only materials with extremely high thermal conduc- 
tivity can survive in these severe environments without melting 
or otherwise losing their structural integrity. Copper and dilute 
copper-base alloys offer the highest possible thermal conduc- 
tivity with mechanical properties sufficient to withstand static 
and cyclic loading. 

Mechanical durability of copper chambers can be limited by 
several potential failure modes: creep-rupture, low-cycle ther- 
mal fatigue, and thermal-mechanical ratcheting. The useful life 


of the throat region of the space shuttle main engine, a rectan- 
gular channel chamber, is governed by thermal-mechanical 
ratcheting. Ratcheting takes the form of local bulging and 
thinning of the rectangular channel wall separating the hot-gas 
and coolant sides, as indicated by the “doghouse” configuration 
in figure 1 (from ref. 1). The rate of ratcheting and wall thinning 
is a strong function of the severity of thermal loading (both 
maximum temperature and thermal gradients), the thermal and 
structural compliance, and the geometry of the combustion 
chamber and is a lesser function of the material strength and 
ductility. Elimination of the ratcheting failure mode by im- 
proving material strength and ductility is highly unlikely 
because extremely large increases in material properties are 
required. Only incremental improvements in material proper- 
ties are feasible at this time. However, as seen from the results 
of the current experimental and analytical study, significant 
improvements can be made in the geometry of the cooling 
channels, and these considerably reduce the amount of ratcheting 
and thereby extend life. Unfortunately, the type of failure ex- 
pected, thermal-mechanical ratcheting, cannot be analytically 
modeled within the present state of the art. Both thermal- 
mechanical ratcheting and fatigue are strainrange dominated; 
therefore, it is reasonable to infer that a fatigue analysis would 
provide a fair estimate of the relative lives of channel geometries. 

In this study the structural durabilities of a tubular channel 
chamber (fig. 2) and a rectangular channel chamber were 
compared. Experimentally, subscale hardware was designed 
and tested to isolate the effects of geometry on thrust chamber 
life. The experimental hardware was analytically modeled. A 
new viscoplastic, constitutive stress-strain model for copper 
was incorporated into the chamber models, and the two con- 
figurations were compared at nominal operating conditions. 
The experimental results from subscale chamber testing con- 
firmed the expectations of the analyses, that significant im- 
provements in cyclic durability can be expected from more 
compliant, less constrained channel constructions. 
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Figure 1 . — Milled rectangular channels, (a) Before repeated 
firings, (b) After repeated firings. 



Figure 2. — "Booked" tubular channel construction. 


Experimental Investigation 

The testing of both chamber designs has been previously 
reported. The rectangular channel testing is reported in refer- 
ence 1, and the tubular channel testing is reported in refer- 
ence 2. The apparatus, test procedure, and results are summarized 
here. 


Apparatus 

A broad base of data on copper-alloy combustion chambers 
with rectangular channel coolant passages has been obtained on 
a low-cost, subscale rocket engine test apparatus at NASA 
Lewis Research Center. The apparatus is used to study the 
cyclic durability of rocket chamber walls, to screen candidate 
chamber-wall materials, and to evaluate fabrication techniques 
for improved life and performance. The test apparatus, shown 
in figure 3, consists of an annular injector; a 15.24-cm (6.0-in.) 
long, liquid-hydrogen-cooled outer cylinder, which serves as 
the test chamber; and a water-cooled centerbody, which forms 
the combustion, sonic throat, and expansion sections of an 
annular rocket engine. It is at the sonic throat station that the 
maximum heat flux is encountered, and this is the location that 
is studied for fatigue. The standard test chamber configuration 
has 72 rectangular coolant passages (ref. 1). 

A variation of this standard configuration was used for the 
comparison reported herein. A new test chamber was designed 
with a stacked bundle of 72 oxygen-free, high-conductivity 
(OFHC) copper tubes as the coolant passages, as shown in 
figure 4. The same number of coolant passages is used in the 
standard rectangular channel chambers, thus allowing direct 
comparison of results. To provide a coolant flow area equiva- 
lent to the rectangular channel configuration and to maintain 
the same wall thickness, round tubes were required to be 
“booked” (i.e., compressed to form the oval shape shown in 
figure 5). 

Procedure 

The test procedure for the tubular channel chambers was 
identical to that used for fatigue tests of the rectangular channel 
chambers described in reference 1. Gaseous hydrogen and 
liquid oxygen were used as propellants at a nominal mixture 
ratio of 6.0 at a chamber pressure of 4 1 37 kPa (600 psia). The 
maximum hot-gas-side wall temperature was controlled by a 
separate mass flow of liquid hydrogen through the coolant 
passages. The chamber was fired for 1 .7 sec to allow it to reach 
thermal equilibrium and then shut down for 1.8 sec while the 
coolant continued to flow to achieve steady-state cold condi- 
tions. The total cycle time was 3.5 sec, as shown in figure 6. 
This cycle was repeated approximately 50 times. The combus- 
tion chamber was then inspected and subjected to another 50 
cycles. This sequence was repeated until a combustion cham- 
ber failure was detected by sensing a coolant passage leak. 

Results 

Two tubular channel chambers were tested for comparison 
with the rectangular channel chambers of reference 1 . The first 
chamber was test fired at the nominal combustion conditions 
and a hot-gas-side wall temperature of 570 °C (1060 °F). 
Failure occurred after 637 cycles. The second chamber was 
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Figure 4. — Stack of OFHC copper tubes and completed 
subscale thrust chamber. 


also test fired at nominal combustion conditions but was run 
with a hot-gas-side wall temperature of 671 °C ( 1 240 °F). This 
chamber failed after 256 cycles (ref. 2). 

Rectangular channel chambers were tested at the same 
nominal operating conditions as the tubular chambers and at 
a wide range of hot-gas-side wall temperatures (ref. 1). The 



Figure 5. — Diagram of tube dimensions (round to booked). 


number of cycles to failure expected for a rectangular channel 
chamber of this configuration tested at 538 °C (1000 °F) is 
250 cycles and 1 17 cycles for one tested at 671 °C (1240 °F) 
(ref. 1). Comparing the test fire results for the two types of 
chambers shows two to three times more cycles to failure for the 
tubular chamber. 
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Figure 6. — Schematic cycle of loading (typical cycle history). 

Analytical Investigation 

Finite Element Thermal and Structural Analyses 

Accurate predictions of the cyclic lives of the tubular and 
rectangular channels require accurate assessments of the in- 
elastic response of these channels. Advanced viscoplastic 
models treat all the inelastic strains (including creep, plastic, 
relaxation, etc.) as a single entity and thus automatically in- 
clude any interactions among them. Viscoplastic models, 
therefore, provide realistic descriptions of inelastic response. 
In this study a viscoplastic model put forth by Freed and Verrilli 
(ref. 3) was employed to calculate the inelastic response of the 
channels. This viscoplastic model is described briefly in the 
appendix. The material constants for copper for this model are 
given in table I. 

The constitutive equations of the viscoplastic model are 
highly nonlinear and mathematically “stiff.” Applying visco- 
plastic models to structural engineering problems, such as the 
present channel comparison problem, requires the develop- 
ment of suitable numerical solution technologies employing, 
for example, the finite element method. In the present work a 
numerical solution technology put forth by Arya and Kaufman 
(ref. 4) was employed. The general-purpose finite element 
program MARC (ref. 5) was used to perform the finite element 
calculations. Time integration of the stiff constitutive equations 


TABLE I.— MATERIAL CONSTANTS FOR COPPER 
[From reference 3.] 


Elastic modulus, £, MPa 

165 000- 1257, Tin °C 

Coefficient of thermal expansion, a, °C -1 . 

... 15x10- 6 + 5x10- 9 T, T in °C 

Poisson’s ratio, v 

0.34 

Constant, A, s -1 

50x1 0 6 

Constant, D 0 , MPa 

1.5 

Constant, h, MPa 

500 

Constant, //, MPa 

5000 

Limiting value of backstress at kinematic 
saturation, L, MPa 

25 exp(-77300) 

Constant, n 

5 

Activation energy, Q, J/mole 

200 000 

Constant, 5, MPa 

14.3 

Melting point, T m , K 

1356 


532 Elements, 

601 nodes 

Figure 7. — Finite element models of 
rectangular and tubular channels. 

of the viscoplastic model requires efficient and accurate inte- 
gration techniques. Because the use of an explicit integration 
method does not require the assembly and inversion of Jacobian 
matrices, a time integration strategy, developed by Arya et al. 
(ref. 6) and based on the explicit forward Euler method, was 
employed here. This integration strategy has been successfully 
applied to numerous problems and has been found to work well 
for problems that involve complex thermal and/or mechanical 
loadings (refs. 7 and 8). 

The finite element models for the rectangular and tubular 
channels were generated by using the PATRAN program 
(ref. 9). These models are shown in figure 7. Generalized plane 
strain, isoparametric elements were used to model the smallest 
repeating segment of the chamber walls. The use of these 
elements facilitated quasi-three-dimensional analyses for the 
channels. The rectangular channel model consisted of 532 
elements and 601 nodes, whereas the tubular channel model 
was constructed with 552 elements and 6 1 5 nodes. These finite 
element models were used to perform the thermal and structural 
analyses presented here. 

Steady-state heat transfer analyses were performed to obtain 
the temperature profiles for the rectangular and tubular chan- 
nels. Relevant properties and boundary conditions were input 
to the MARC program. The steady-state temperature distribu- 
tions for the rectangular and tubular channels obtained from 
these analyses, shown in figure 8, were akin to each other. 

The cyclic thermal and mechanical loadings used in the finite 
element analyses are shown in figure 6 and table II (also see 



552 Elements, 
615 nodes 
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After first cycle 



After fifth cycle 



Figure 8. — Temperature distribution in rectangular and 
tubular channels during fifth cycle. 


Figure 10. — Effective stress in tubular channel after first 
and fifth cycles. 
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Figure 9. — Effective stress in rectangular channel after 
first and fifth cycles. 


Figure 1 1 . — Maximum effective strain in rectangular and 
tubular channels during fifth cycle. 







TABLE II.— CYCLIC TEMPERATURE 
AND PRESSURE LOADING HISTORIES 
[From reference 10.] 



Wall pressure, kPa (psia) 

Cold phase 

Hot phase 

Hot-gas side 
Coolant side 

96.5(14.0) 

5100(740) 

2780 (403) 
6550 (950) 


ref. 10). The temperatures during the transient stage of the 
thermal cycle were obtained from corresponding steady-state 
temperatures by employing linear interpolation. The duration 
of a typical thermal-mechanical cycle was 3.5 sec. Time- 
varying nodal temperatures and elemental pressure loadings 
were applied to trace the loading cycle shown in figure 6. The 
variations in material properties with temperature (see table I) 
were taken into consideration for these computations. 

The effective stress 1 distributions in the rectangular and 
tubular channels are shown in figures 9 and 10, respectively. 
These figures exhibit effective stress in the channels at the end 
of the first and fifth loading cycles. Examination of these 
figures reveals a significant redistribution of effective stress 
due to inelastic deformation. In the case of the rectangular chan- 
nel, for example, at the end of the first cycle the effective stress 
in the lower part of the segment was higher than that in the upper 
part. But at the end of the fifth cycle, the effective stress in the 
lower part of the segment was lower than that in the upper part. 
Similar conclusions can be drawn for the tubular channels. 

Figure 1 1 shows the maximum effective strains 1 in the 
rectangular and tubular channels during the fifth loading cycle. 
The maximum effective strain in the tubular channel was 
considerably lower than that in the rectangular channel. Note 
also from this figure that the location of maximum effective 
strains differed for the rectangular and tubular channels. 

The deformed shapes of the rectangular and tubular channels 
are plotted in figure 12 at the end of the fifth cycle. The 
deformations in the two channels are magnified by a factor of 
1000 to facilitate visual interpretation. The figure shows a 
significantly larger deformation of the rectangular channel than 
of the tubular channel. 

Figures 1 3 and 1 4 display the thinning of channel walls in the 
circumferential (X) and radial ( Y) directions. The values plotted 
are at the ends of five loading cycles. The figures reveal that the 
rectangular channel wall thinned at a more rapid rate than did 
the tubular channel wall. For example, at the end of the fifth 
cycle the thickness of the rectangular channel wall had been 
reduced approximately twice as much as the tubular channel 


1 Effective stress : 





Effective strain : 




Figure 1 2. — Deformed shapes of rectan- 
gular and tubular channels after five 
cycles (magnification factor = 1 000). 



Cycles 


Figure 13. — Thinning in X-direction of rectangular and 
tubular channel walls after different numbers of cycles. 


wall thickness. Looking at the deformed shapes in figure 1 2 the 
thinning in the radial direction is quite obvious. 

It is clear from this discussion that for the same thermal- 
mechanical loading cycles: (1) the deformation of the tubular 
channel was strikingly smaller than that of the rectangular 
channel and (2) the thinning of the tubular channel wall in the 
circumferential and radial directions was much slower than that 
of the rectangular channel wall. It may thus be concluded that 
using a tubular geometry for the coolant channel can greatly 
increase the durability and life of thrust chambers. 
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Cycles 

Figure 14. — Thinning in /-direction of rectangular and 
tubular channel walls after different numbers of cycles. 

Life Analyses 

Life analyses of the rectangular and tubular channel com- 
bustion chambers were performed in accordance with an estab- 
lished procedure referred to as the “local strain” approach. The 
principal elements of this approach are 

(1) Thermal and mechanical (pressure) mission loads imposed 
on the chamber 

(2) Heat transfer coefficients for the fluids and the structural 
components 

(3) Thermal conductivity and thermal expansion properties of 
the material 

(4) Cyclic elastoplastic constitutive stress-strain-temperature- 
time behavior of the material 

(5) Low-cycle fatigue behavior of the material as represented 
by the cyclic effective plastic strainrange versus the cyclic 
crack initiation life 

The first four elements are used as input to finite element 
analyses, which produce an output of local stresses, strains, and 
temperatures for each element as a function of mission opera- 
tion time. Scanning the output locates the critical finite ele- 
ments where temperature extremes and strainrange are largest 
and which hence are the most likely locations for thermal low- 
cycle fatigue cracks to initiate. 

The material in the critical locale is deemed to have failed 
from fatigue crack initiation when a crack of critical size 
develops. The critical crack size is one large enough to cause 
fatigue separation of a uniaxial-loaded laboratory coupon sub- 
jected to an equivalent level of cyclic strain, temperature, and 
exposure time. No account is taken herein of cyclic propaga- 
tion of cracks larger than this somewhat arbitrary critical 
initiation size. Low-cycle, strain-life fatigue curves for the 
copper material were approximated from data available for 
related thermal and temporal variations in cyclic conditions. 

The thermal-mechanical loading cycle shown in figure 6 was 
applied repeatedly to both chamber designs. The maximum 


TABLE III.— TEMPERATURES AND RANGES OF STRESS AND 
STRAIN AT CRITICAL LOCATIONS IN RECTANGULAR 
AND TUBULAR CHANNELS 


Channel 

Metal temperature, 
°C (°F) 

Effective stress 
range, 

MPa (psi) 

Effective 

strainrange, 

percent 

Max. 

Min. 

Total 

Plastic 

Rectangular 

Tubular 

538(1000) 
445 (833) 

-218 (-360) 
-218 (-360) 

167 (24 221) 
210(30 458) 

2.60 

1.05 

2.46 

.87 


and minimum temperatures, effective strains, and effective 
stresses at the most critical location for each chamber design are 
presented in table III. This information evolved from the 
thermal and structural analyses of the chambers as they were 
cycled between the temperature extremes of the mission load- 
ing cycle. Supplemental calculations indicated that the pres- 
sure loading, acting alone, would have a negligible effect on the 
deformation and life of the chambers, (i.e., the pressure- 
induced stresses were well below the yield strength of the 
material at the nominal maximum operating temperature of 
538 °C (1000 °F)). The associated creep rates also were small. 
Despite the small pressure influence when acting alone, the 
simultaneous application of the pressure loading and the ther- 
mal gradients gave rise to large cyclic strains and an accelerated 
rate of ratcheting. The thermal gradient (and associated ther- 
mal constraint) contribution to the range of total strain at the 
critical location was 2.6 percent for the rectangular channel 
design (of which 2.46 percentage points were plastic strain) and 
1.05 percent for the tubular channel design (of which 0.87 
percentage point was plastic strain). The time per cycle was 
extremely short (3.5 sec), and only 2.3 sec were spent at 
temperatures above the minimum cryogenic level. Conse- 
quently, during heating and cooling, the average strain rate was 
quite high (on the order of 1 to 2 percent per second). Because 
of the high strain rate, insufficient time was available for 
appreciable time-dependent, thermally activated creep defor- 
mation to occur. Any creep that might be present would occur 
most likely as stress relaxation during the brief, 0.8 sec, hold 
period at maximum temperature. Local stress relaxation of as 
much as 10 percent from the peak compressive stress would 
introduce a compressive creep strain of less than about 0.01 
percent. This strain is approximately three orders of magnitude 
smaller than the plastic strainrange and would not be expected 
to affect the cyclic life of the channel wall. Consequently, for 
nominal chamber operation there was little concern for creep- 
fatigue interaction during hot-fire testing of the chambers, and 
the life prediction analysis was simplified to a simple cyclic 
elastoplastic condition. 

Pertinent high-temperature, low-cycle, plastic strain fatigue 
data have been reported by Conway et al. (ref. 1 1 ) for annealed 
OFHC copper at 538 °C (1000 °F) in a nonoxidizing argon 
atmosphere. These data are shown in figure 1 5 as total effective 
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Figure 15. — Low-cycle fatigue resistance of OFHC copper 
at 538 °C (1000 °F) in argon at cyclic strain rate of 0.2 
percent per second. 


plastic strainrange versus cycles to crack initiation. The data 
had been generated at a cyclic strain rate of 0.2 percent per 
second, which is about one-tenth the rates experienced in 
the chambers. At this cyclic strain rate the extent of time- 
dependent, thermally activated cyclic creep strain is negligible. 
Hence, the data meet an important criterion for applicability to 
the current analyses — negligible creep-fatigue interaction. 
Another important criterion for using Conway’s isothermal 
data is their applicability to the variable temperatures experi- 
enced throughout the mission cycle of the combustion cham- 
bers. Because the tensile ductility of OFHC cooper is virtually 
insensitive to temperature, it is expected (ref. 12) that the low- 
cycle fatigue resistance (plastic strainrange versus life) would 
also be insensitive to temperature. Tensile test reduction of area 
varies from 81 to 98 percent over the temperature range of 
concern (approximately -218 to 538 °C; -360 to 1000 °F). If 
the fatigue failure relation is independent of temperature, and 
only one deformation mechanism is present (cyclic plasticity), 
it is reasonable to assume that the fatigue failure relation will be 
directly applicable to the thermal cycling conditions being 
considered. In other words, one curve, determined at any 
isothermal temperature within the cyclic range of -218 to 
538 °C (-360 to 1000 °F) should be applicable to all tempera- 
tures within the range if the strains in the cycle are predomi- 
nately plastic. The latter was the case for both the isothermal 
fatigue test data and the strains occurring at the critical locations 
within the chambers. Note that environmental oxidation inter- 
action was precluded by use of an inert argon atmosphere for 
coupon testing and by the hydrogen-rich, reducing atmosphere 
present during chamber testing. Note also that OFHC copper 
is not embrittled by the presence of hydrogen (ref. 13). 

For temperatures far above 538 °C ( 1 000 °F), however, creep 
and associated creep damage cannot be ignored because they 
increase with temperature at an exponential rate. Some of the 
chamber tests experienced temperatures as high as 671 °C 
(1240 °F); therefore, creep strain due to the internal channel 
pressure acting alone would become (based on a simple Larson- 
Miller parameter calculation) nearly three orders of magnitude 


greater than at 538 °C (1000 °F). Consequently, the time- 
dependent deformation would become a nonnegligible fraction 
of the plastic strain with a resultant loss in cyclic life expect- 
ancy. Quantitatively acceptable calculations of the creep-strain 
contributions were not possible without performing three- 
dimensional finite element viscoplasticity analyses for the 
highest temperature chamber tests. Qualitatively, a large 
decrease (approaching an order of magnitude in life) in cyclic 
durability could be expected. 

Predicting fatigue crack initiation life for the tubular and 
rectangular chamber designs under nominal thermal and me- 
chanical loading conditions is a rather simple task. The two 
equivalent plastic strainranges of 2.46 and 0.87 percent were 
entered into the Manson-Coffin law shown in figure 1 5, and the 
respective cyclic lifes of 70 to 750 cycles to failure were 
determined. The greater fatigue life of the tubular relative to the 
rectangular channel design is significant. A more than tenfold 
improvement in predicted low-cycle fatigue life is due directly 
to the reduction of the thermally induced strainrange in the 
tubular channel design. 


Discussion 

This effort showed that the strain levels for the tubular 
channel were significantly lower than those for the rectangular 
channel. Strainrange was used to predict a significant increase 
in mission life. The analyses predicted a reduction in strainrange 
from 2.6 percent in the rectangular channel to 1 .05 percent in 
the tubular channel. Such a dramatic reduction in strainrange 
was expected to provide a tenfold increase in life by “the local 
strain approach.” The test results, however, only showed a two- 
to threefold increase in life. 

Both analyses and tests showed that significant improve- 
ments in life can be achieved by accommodating strain instead 
of resisting it, but the amounts of improvement differed. This 
difference was principally the result of our inability to accu- 
rately include all effects of the real environment. Because of 
this reality, only relative comparisons could be made. The 
difference in life indicated by the analyses and the tests was due 
primarily to two reasons. 

First, a fatigue analysis was performed to estimate the lives 
of the tubular and rectangular channel designs. Thermal- 
mechanical ratcheting, not fatigue, is the failure mode of these 
types of chambers. But both fatigue and thermal-mechanical 
ratcheting are dominated by the magnitude of the strainrange. 
Therefore, it is reasonable to assume that a fatigue analysis 
would provide a good, but not exact, indicator of relative life 
improvement. 

Second, there was a difference in operating versus assumed 
temperatures. The analyses used a maximum wall temperature 
of 538 °C (1000 °F) for both chamber configurations so that a 
direct comparison could be made. In testing the two configu- 
rations it was difficult to achieve an exact temperature of 
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538 °C (1000 °F). The tubular channel was run slightly hotter 
than the analysis temperature (on the order of 571 °C(1060 °F)). 
Consequently, some of the potential improvement in durability 
because of a more compliant geometry was offset by a higher 
temperature with attendant higher rates of damage due to the 
probable introduction of cyclic creep strain. Furthermore, 
because of the nominally higher temperature of the tubular 
channel, greater surface roughening was observed. This phe- 
nomenon, known as blanching (ref. 14), escalated the tubular 
wall temperature to an estimated value of 830 °C (1525 °F) at 
the time testing was suspended. This additional increase in 
temperature of the tubular channel chamber reduced its life 
even further relative to both the tubular channel analysis and the 
experimental results of the rectangular channel chamber. 

It is obvious from the results of both the testing and analyses 
that a tubular channel geometry reduces the strainrange and 
therefore increases the life of a rocket thrust chamber. 


Conclusions 

Presently, advanced hydrogen-oxygen rocket combustion 
chambers have a limited life due to thermal-mechanical fatigue 
and ratcheting. This ratcheting is a strong function of the 
material strength and the strainrange. In the past numerous 
attempts have been made to extend the useful life of rocket 
combustion chambers by improving the chamber material 
strength. These efforts resulted in little life improvement 
because the strain levels at the high heat fluxes present in rocket 
combustion chambers are in the range 2 to 4 percent. This level 
of strain is roughly an order of magnitude greater than the yield 
strain of copper or copper alloys. 

The work herein takes steps to dramatically improve com- 
bustion chamber life by reducing the strainrange the material 
experiences. This reduction is accomplished by modifying the 
geometry of the hot-gas-side surface of the combustion cham- 
ber. The tubular geometry reduces the compressive strain, thus 
reducing the total strainrange and extending life. 

The viscoplastic analysis, durability analysis, and test results 
reported herein verified that the more compliant tubular con- 
struction extends the life of the combustion chamber. The 
viscoplastic analysis found that the total strainrange of the 
tubular channel was reduced to 1.05 percent from 2.6 percent 
in the rectangular channel. The durability analysis found that 
the reduction in strainrange could extend the life of the combus- 
tion chamber by as much as 1000 percent. The test results found 
that the compliant chamber extends life by 300 percent. 


D. The kinematic, or flow-induced, anisotropic hardening is 
accounted for by the backstress B y, and the isotropic hardening 
is incorporated through the drag stress D. There is no yield 
surface in the model. The assumption of uncoupled dynamic 
and static recovery terms provides a simpler mathematical 
structure to the model. 

The total strain rate £y is assumed to be the sum of elastic S- , 
inelastic ejj 1 , and thermal strain rate components. This can 
be written as 

*«/=$ + #+*? (U = 1.2,3) 

Here e denotes strain and a dot over a symbol indicates its 
derivative with respect to time t. Note that the inelastic strain 
in a viscoplastic formulation includes plasticity, creep, relax- 
ation, etc., and their interactions. 

The material is assumed to be isotropic, and Hooke’s law is 
used to relate the stress rate Ofy to the elastic strain rate sfj , that 
is, 


• el 1 + V . V . - 

ij p ° ij p G kk ij 


In the preceding equation, E is Young’s modulus, v is 
Poisson’s ratio, and 8^ is the Kronecker delta. Following the 
usual summation convention, the repeated subscripts in the 
preceding equation and elsewhere imply summation. 

We define deviatoric stress Sy as 


S ij - °ij 3 a kk S ij 


Denoting the backstress as By, the effective stress X ; y is defined 
as 


^ij - $ij &ij 

The second variant of effective stress J 2 is 


and the function F is defined by 


F = 


_ 2 _ 

D 


where D is the drag stress. 

The constitutive equations of the model now follow. We first 
define the thermal diffusivity function G, which contains the 
temperature dependence of the model. It is given by 


Appendix — Viscoplastic Model 

The viscoplastic model put forth by Freed and Verrilli 
(ref. 3) and used in the present work incorporates two internal 
state variables — tensorial backstress By and a scalar drag stress 
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in which K is the Boltzmann constant, Q is the activation 
energy, T is the absolute temperature, and T m is the melting 
point of the material. 

The function Z, called the Zener-Hollomon parameter, has 
the following expressions: 

Jaf" f<i 

^ | A exp[n(F - 1)] F > 1 

The flow law is now given by 


The function r is constrained by the following dissipativity 
condition: 


r > Z 
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r b\ yi 
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Fh — — 
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ID 


v 'J 


where 


A 


— B B 
2 lJ lJ 



The differential equations governing the growth of the inter- 
nal state variables B and D are 


B u= H 


e ,n 

IJ 



\ 

J 


and 


where 


D = h 



I 2 =0Z 

In the above equations H , L, and h are inelastic material 
constants. The constant L denotes the limiting value of the 
backstress at kinematic saturation (ref. 3). 

The following expressions then define the recovery func- 
tion r 


r(G) = 



D = D 


O 

D> D 


Jag" -1 g<i 

I A exp[n(G- 1)]/G G>1 


and 


G = 


L 

S-D 


where S and D () are material constants. 


This condition must be satisfied for the theory to be thermody- 
namically admissible. The material constants for this model 
characterized for copper are listed in table I. 
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